Kondo-type transport through an interacting quantum dot coupled to ferromagnetic 

leads 
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We investigate the equilibrium and out-of-equilibrium Kondo effects in a single-level interacting 
quantum dot connected to two ferromagnetic leads. Within the non-crossing approximation, we 
calculate the total density of states (DOS), the linear conductance, and the nonlinear differential 
conductance for both the parallel and the anti-parallel alignments of the spin polarization orientation 
in the leads, followed by a brief discussion regarding the validity of this approach. Numerical 
calculations show that for the anti-parallel alignment, a single Kondo peak always appears in the 
equilibrium DOS, resulting in the conventional temperature behavior in the linear conductance 
and the zero-bias maximum in the differential conductance. The strength of the DOS peak is 
gradually suppressed with increasing polarization, due to the fact that formation of the Kondo- 
correlated state is more difficult in the presence of higher polarization. On the contrary, for the 
parallel configuration the Kondo peak in the DOS descends precipitately and splits into two peaks 
to form a very steep valley between them. This splitting contributes to the appearance of a "hump" 
in the temperature-dependent linear conductance and a nonzero-bias maximum in the differential 
conductance. Moreover, application of a bias voltage can split each Kondo peak into two in the 
nonequilibrium DOS for both configurations. Finally we point out that the tunnel magnetoresistance 
could be an effective tool to demonstrate the different Kondo effects in different spin configurations 
found here. 

PACS numbers: 72.10.Fk, 72.15.Qm, 72.25.Dc,73.63.Kv 



I. INTRODUCTION 



The discovery of the Kondo effect in a quantum dot 
(QD) connected to two normal reservoirs has stimulated 
much experimental and theoretical interests in this many- 
body phenomenon,"'^ resulting in several novel findings, 
such as splitting of the Kondo peak under nonequilibrium 
condition;^ an unususal enhancement of conductance in 
the cases when an even number of electrons resides in a 
QDj^ and double peaks of the differential conductance 
in coupled double QDsiS^A Recently, an increasing atten- 
tion has been paid, because of its potential application in 
the megnetoelectronics and quantum computer^ to the 
spin-polarized electron tunneling through systems con- 
sisting of two ferromagnetic (FM) leads sandwiched by 
a QD, which features a significant Kondo effect at low 
temperatures when it connects to normal leads. Gen- 
erally speaking, the Kondo resonance in the density of 
states (DOS) at the Fermi energy originates from screen- 
ing of the dot spin due to the exchange coupling with the 
conduction electrons. Therefore it is interesting to ob- 
serve whether the Kondo-correlated state can form when 
the conduction band is of spin polarization, and if so, 
what is the difference from the conventional patterns. 

In a recent paper, Sergueev et alP and Zhang et alSi 
presented a theoretical analysis of the transport char- 
acteristics of such a FM/QD/FM system, using the 
ansatz proposed by Ng^" and the standard equation-of- 
motion (EOM) technique for the retarded Green's func- 
tion. They found that there is always a sharp single 
Kondo resonant peak in nonlinear differential conductiv- 



ity at zero bias, regardless of the polarization orienta- 
tion of the two leads, parallel (P) or anti-parallel (AP) 
configurations. On the contrary, Martinek et ali^ re- 
ported a markedly different result on similar systems, 
based also on an EOM technique and an additional as- 
sumption to replace the bare level of QD in the resultant 
self-energy expression with the one self-consistenly deter- 
mined. They found that for P alignment of the lead mag- 
netizations, the Kondo resonances in the DOS split for 
spin-up and -down electrons, thus the differential conduc- 
tance exhibits nonzero-bias maximum and the linear con- 
ductance drops to a low value even the polarization being 
as small as 0.2. Later on, Lii and LiuiS. reported also the 
similar splitting by applying Ng's ansatz. Ma et alM em- 
ployed the tinite-U slave-boson mean-field approachi^ti^ 
to investigate the spin-polarized transport of this system 
and found drastically different behaviors depending on 
the polarization alignment of the two leads. In view of 
the ongoing controversy, further analysis on the influence 
of the spin polarization on the Kondo-correlated state in 
a QD is desirable, preferably based on more advanced 
schemes. 

In this paper, we employ the non-crossing approxima- 
tion (NCA) based on the auxiliary boson technique to 
carry out a detailed analysis of the FM/QD/FM Kondo 
problem. The NCA is a diagrammatic technique to 
sum all the non-crossing diagrams in the leading order 
\V\'^/N {V is the hopping matrix element between the 
local electron and conduction electrons, and N is the 
number of spin degeneracy of the local level) ^ And it is 
proved to be an accurate approach for the case of iV = 2 
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(of interest for a QD), even under out-of-equalibrium 
conditions For a symmetric FM/QD/FM system 
considered in the present paper, although the individ- 
ual QD is exactly spin degenerate, the spin-polarized 
leads can lift this degeneracy in the case of P alignment 
through the tunneling induced spin-related self-energies. 
However, if the spin polarization strength p of the two 
leads is small enough, the self-energies are of weak spin- 
dependence, and as a consequence the NCA should be 
adequate to treat this systems. Actually, we calculate 
the equilibrium and nonequilibrium transmissions for the 
FM/QD/FM system with the polarization strength as 
low ss, p — 0.2 for the P configuration. Numerical results 
show a complete splitting for the equilibrium transmis- 
sion (DOS) and a nonzero-bias maximum for the voltage- 
dependent differential conductance, which are in qualita- 
tive agreement with those of Martinek-'^^ and Liii^. 

The rest of the paper is organized as follows. In sec- 
tion II we describe the model Hamiltonian and give a 
brief description of applications of the infinite-f/ NCA to 
the spin-dependent transport in FM/QD/FM systems, as 
well as a qualitative discussion about its validity. In sec- 
tion III numerical computations and discussions are pre- 
sented, including the equilibrium and out-of-equilibrium 
DOS, the linear conductance, and the nonlinear differen- 
tial conductance, for both the P and the AP polarization 
configurations. Finally, all the results are summarized in 
Section IV. 



II. MODEL AND FORMULATION 

The system Hamiltonian for a quantum dot with a 
single spin-degenerate energy level {N — 2) connected 
to two ferromagnetic leads is written as 

ri,k.a <y 

+t/nd|nd| + ^ (Vc^fegC^a +H-C.), (1) 

r\.k^<y 

where cj^^^ ic^ka) and cj^^ [cdcj) are the creation (anni- 
hilation) operators for electrons with momentum fc and 
spin a in the lead r] (= L, R) and for a spin-a electron 
on the QD, respectively. The third term describes the 
Coulomb interaction among electrons on the QD, which 
is assumed to be infinite ([/ — > oo) in the present paper, 
forbidding double occupancy. The fourth term represents 
the tunneling coupling between the QD and the reservoirs 
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2<5(co - e^feg), (?7 = L,R). (2) 



In the wide band limit, F^ is assumed to be constant. For 
the identical leads and symmetric barriers, of interest 
in the present investigation, the ferromagnetism of the 
leads can be accounted for by the polarization-dependent 



couplings Ff = Ff = (1 -hp)Fo, F[ = Ff = (1 - p)Fo 
for the P alignment, while Fj' = Ff = (1 + p)Fo, F[ = 
= (1 -_p)Fo for the AP alignment. Fq, and p (0 < 
p < 1) describe the tunneling coupling between the QD 
and the nonmagnetic leads, and the polarization strength 
of the leads. Under this approximation, the current / 
through the QD can be expressed in terms of the total 
transmission T(lo) as^li 



dc^[/i(c.)-/fl(c.)]r(c^). 



where 



T(u;) 



F^F 



^ F£ + F« 



PrT(w), 



(3) 



(4) 



with Paiyj) — ^:^IniG'JJ.(ijj) being the DOS for spin-cr elec- 
trons. GJ^(a;) is Fourier transform of the retarded Green's 
function. 



(5) 



The main purpose of this work is to calculate the DOS 
Pa-(w) as a function of temperature T, bare-level energy 
td and bias voltage V for different polarization config- 
urations and strengths p, and the associated linear and 
nonlinear conductance. 

According to the infinite-C/ slave-boson approach, the 
ordinary electron operators on the QD can be decom- 
posed into a boson operator h and a pseudo-fermion op- 
erator /o-, 



(6) 



with a constraint for the auxiliary operators b^h + 
J2cr fa-fa = 1- thc slavc boson representation, the 
Hamiltonian |^ for the FM/QD/FM systems becomes 



^ ^ ^ykaclka^nka + fafg- 

rj^k.a a 

+ E(v4>/-+H.c.). 

?7 , A: , cr 



(7) 



In order to evaluate the DOS Pa{'-^): Wingreen et alm^ 
generalized the NCA to study the nonequilibrium prop- 
erties of the Anderson model connected with two nor- 
mal conduction bands, using the Keldysh nonequilibrium 
Green's function formalism. It is well-known that the 
NCA is a self-consistent conserving perturbation expan- 
sion for the pseudo-fermion and slave-boson self-energies 
to first order in the effective coupling J = At the 

lowest order in perturbation diagrams the boson self- 
energy involves the bare fermion propagator while the 
fermion self-energy involves the bare boson propagator. 
By replacing these bare propagators with the dressed 
auxiliary particle propagators in the Feynman diagram, 
one can obtain a set of coupled integral equations, which 
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self-consistently determine the self-energies of these aux- 
ihary particles. Solving these coupled equations is equiv- 
alent to summing up a subset of diagrams to all orders in 
J. Furthermore, it can be proved that the NCA includes 
all diagrams of leading orders in 1/A^,^^ Therefore, the 
NCA is expected to be a quantitative approach in the 
limit of large N. For a QD connected with normal leads, 
A'^ = 2, it is already proved to be satisfactory in quali- 
tatively describing the linear and nonlinear Kondo-type 
transport 

Unfortunately, when the level degeneracy is broken, 
the NCA could produce spurious peaks in the DOS, thus 
is unreliable for transport investigation. For example, 
as mentioned by Wingreen/* the NCA without vertex 
corrections produces an additional Kondo peak at the 
chemical potential in a finite magnetic field due to a 
false sefl-interaction of each level, whereas other meth- 
ods find that the Kondo peak splits into two peaks. Sim- 
ilarly, the NCA without vertex corrections seems to be 
also inappropriate for the FM/QD/FM systems, because 
spin-related tunneling lifts the level degeneracy in the 
QD. However, the present situation is somewhat differ- 
ent from magneto-transport. The degeneracy lifting is 
evident in the presence of a magnetic field, while is de- 
pendent on the relative polarization orientation of the 
two leads and of course the strength of polarization for 
the FM/QD/FM systems. 

When the polarization orientations of the two FM leads 
is anti-parallel, the self-energies are actually independent 
of spin and the degeneracy remains as two for the case 
of the identical leads and symmetric barriers. On the 
other hand, for the P configuration spin-related tunnel- 
ing results in the self-energies to be different for spin-up 
and -down electrons. However, it is natural that the de- 
viation depends on the polarization strength p. Namely, 
the NCA could still be rehable for the FM/QD/FM sys- 
tems with small enough p. Numerical calculations in the 
next section show that the conventional Kondo peak in 
DOS indeed splits into two peaks completely at the po- 
larization p — 0.2 for the P configuration. Even though 
one can observe an additional peak located at the chem- 
ical potential but an order of magnitude smaller than 
the two real Kondo peaks, these results are in agreement 
with previous predictionsiii Consequently, this convinces 
us that the NCA provides an appropriate description for 
the symmetric FM/QD/FM systems in either the AP 
alignment with arbitrary polarizations or the P configu- 
ration with weak polarizations. Of course, this scheme's 
validity being dependent on the polarization p should be 
carefully checked, for example, by the modified NCA in- 
cluding vertex corrections. This, however, entails numer- 
ically a much heavier task than the original NCA and is 
beyond the purpose of the present paper. We leave this 
examination to a future publication. 

We outline the formulation employed in this paper as 
follows. The interested reader can refer to Ref 18 and 
ReflTgl for detail. In the slave boson representation, the 
retarded Green's functions for the boson and pseudo- 



fermions are defined as 



1 



(8) 
(9) 



with the corresponding retarded self-energies H''(a;) and 
S^^(a;). Furthermore, the "lesser" Green's functions for 
the boson and fermions are related with the "lesser" self- 
energies n<(a;) and T,J^{uj) as 



D<{uj) = D''{lu)U<{uj)D''{lu), 



(10) 

G%{u:) - G%{u)Y.%(u,)G)„(u,). (11) 

The self-consistent NCA equations for out-of-equilibrium 
are 



n''H = E ^ defies uj-f^,)G}Ae), {^'2) 

?7=L,R,(T 

^/-H = E ^ [defies cJ + ^^,)D^ie), (13) 
= E ^ ^defie~u; + ^^,)G%ie),iU) 

where f{x) — [exp(/3a;) + 1]^-^ {(3 — l/k^T) is the Fermi 
distribution function and /i^ is the chemical potential of 
the ry lead. After solving this set of self-consistent equa- 
tions, the imaginary part of the retarded local Green's 
function, the DOS p„{lj), can be calculated within the 
NCA as 

P.H = 1/ ^[D<{e)lu,G)Ae + u) 

+ G<^{e)lmD^{e-u;)], (16) 



where 



Z = j de[D<{e) + Y,G%{e)]. 



(17) 



Finally we can use Eqs.© and to calculate the cur- 
rent through the QD. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Density of State 

In this section we present numerical calculations and 
discussions. First, we deal with the total equilibrium 
and out-of-equilibrium DOS (transmission) Eq.l^J for 
the FM/QD/FM systems with a fixed bare- level energy 
td = —4.0 (Po is used as the energy unit throughout the 
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rest of the paper). It is worth noting that the systems 
considered here belong to the deep Kondo regime and are 
appropriate to demonstrate the strong correlated effects. 

Fig. 1(a) shows the total DOS in the AP configuration 
for several different polarizations p = 0, 0.2, and 0.4 as 
well as various temperatures T = 0.01, 0.02, and 0.04. 
Clearly, a significant Kondo peak remains at the Fermi 
energy (which is chosen to be the the energy zero) un- 
der the addition of the spin-polarized leads. The con- 
tribution of spin-polarized leads is to suppress both the 
Kondo peak and the single particle excitation peak (see 
overall shapes of the DOS in the inset of Fig. 1). This 
suppression is more pronounced with increasing polar- 
ization p. These results are understandable with the aid 
of the following considerations. Suppose that in the ex- 
treme case of two completely spin-polarized leads p = 1, 
spin-down electrons are completely absent in the left lead, 
but electrons in the right lead arc all spin-up and could 
provide compensation to screen the dot spin and to guar- 
antee the formation of the Kondo-corrolatcd singlet state. 
Thus the Kondo peak still exists with a reduced ampli- 
tude. Naturally, the cases of weak polarizations p < 1 
are more likely to form the Kondo state. Note that in- 
creasing temperature can broaden the peak and suppress 
the Kondo resonance as usual. In addition, effects of 
the external bias voltage on the out-of-cquilibrium DOS 
are plotted in Fig. 1(b). For convenience, we choose a 
symmetric voltage drop such that the chemical poten- 
tial /ii = —fJ.R = eV/2 for the left and right leads. As 
expected, we find a splitting of the Kondo peak with a 
width nearly equal to the bias voltage applied between 
the source and drain leads. 

In short, the QD connected with two AP magnetized 
leads develops the same Kondo resonance as the QD with 
two normal leads, whereas the former is suppressed to 
some extent depending on the polarization p. In con- 
trast, the situation is drastically different for the P con- 
figuration as shown in Fig. 2, where we plot the total 
DOS for the polarization p = 0.2 with various temper- 
atures (a) and bias voltages (b). It is clear that the P 
polarization significantly changes the DOS of the QD in 
comparison with the case of p = 0. The Kondo resonance 
splits into two distinct peaks with different amplitudes. 
One moves from the original location w = 0, the Fermi 
energies of the two leads, to a lower energy position, the 
other shifts to the opposite direction at the expense of 
its height. Moreover the magnitudes of both Kondo res- 
onances are largely suppressed by the introduction of spin 
polarization. We observe that there appears a remnant 
peak located at the Fermi energy, which is produced by 
the NCA calculation without vertex corrections as men- 
tioned above. But its amplitude is one order smaller than 
the two shifted Kondo peaks and can be neglected. As a 
result, a very steep valley is found between the two peaks 
with a nearly vanished bottom. Finally we find that in- 
creasing temperature can not only smooth and broaden 
the peaks as usual, but also raise the bottom of the deep 
valley gradually, which can result in a peculiar temper- 



ature dependence of the linear conductance as shown in 
Fig. 3 (in the next subsection). 

Fig. 2(b) depicts effects of changing bias voltage on the 
total DOS (transmission probability) for the same system 
as in Fig. 2(a). If we keep the temperature low, T = 0.01, 
and increase the bias voltage, the two resonances first 
experience suppression and then each of them splits into 
two distinct peaks. Each pair of the peaks has a width 
about equal to the bias voltage. Increasing the temper- 
ature would eventually wash out the peak splitting and 
recover a single but much less pronounced peak in both 
equilibrium and out-of-equilibrium cases. 



B. Linear and Nonlinear Conductance 

In Fig. 3 we plot the calculated linear response con- 
ductance G vs. logT for various polarizations p = 
0,0.1,0.15, and 0.2. In the AP configuration (thin line 
in Fig. 3) the linear conductance Gap exhibits the simi- 
lar overall temperature dependence with those of non- 
magnetic leads p = 0, though smaller in magnitude 
resulted from the suppression of the Kondo resonance 
as addressed in Fig. 1(a). For the opposite orientation, 
the conductance Gp depends strongly on the polariza- 
tion strength p and exhibits no universal T behavior 
at low temperatures. Increasing p, Gp is largely sup- 
pressed first and develops a "hump" as a function of 
temperature. This peak is due to the fact that the 
Kondo resonance is shifted away from the Fermi energy as 
shown in Fig. 2(a). To demonstrate the dramatic change 
of conductance imder different polarization orientations, 
we plot the tunnel magnetoresistance (TMR) defined as 
TMR = {Gp - Gap) /Gap in the inset of Fig. 3. We find 
the TMR arrives at a value as large as 100% at the low- 
est temperature calculated in this work and falls rapidly 
with increasing T, which can be attributed to the pe- 
culiar Kondo resonance in the P configuration. At high 
temperatures Gp approaches the value of nonpolariza- 
tion p = 0, leading to a saturated and small positive 
hnear TMR. 

Fig. 4 shows the linear conductance as a function of 

bare-level energy e^, which can be timed via the external 
gate voltage, for nonpolarization p = 0, as well as the P 
and the AP configurations with p = 0.2. As expected. 
Gap shows the same trend with that of nonpolarization 
though a smaller amplitude. The peak of Gp, however, 
shifts towards the Fermi energy. This means that the 
linear TMR changes its sign at a certain level energy 
and has approximatively a symmetric shape around this 
point. Near the empty orbital regime r,; ~ and near 
the deep Kondo regime ^ —4, the linear TMR reaches 
its maximum value as large as 30% in the temperature 
T = 0.1. 

Nonlinear differential conductance dl/dV is believed 
to be a very useful and sensitive tool in experiments to 
detect the formation of the Kondo-correlated state due to 
its proportionality to T{eV) derived from the current for- 
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mula Eq.Q, assuming that the total transmission T{lu) 
(nonequihbrium DOS) is unchanged under the external 
bias voltage V. So we illustrate in Fig. 5 the calculated 
dl/dV under the P (a) and the AP (b) configurations, as 
well as the nonlinear TMR (c) at various temperatures. 
As pointed out above, because electrons with spin-up and 
spin-down are equally available in the AP configuration, 
the formation of the Kondo-correlated state should not 
be affected. As a result, all the curves in Fig. 5(b) ex- 
hibit a single zero-bias peak and rapid decrease in peak 
height with increasing temperature. When the magneti- 
zation is rotated to the P alignment, adding bias voltage 
can greatly enhance dl/dV at low temperatures. This 
is apparently due to the complete splitting of the Kondo 
peak in the DOS shown in Fig. 2. The nonzero-bias max- 
imum in dl/dV is in good agreement with previous EOM 
calculations fii except for the reduced width of splitting 
and the fine shape in the differential conductance. These 
inconsistences can be attributed to the significant change 
of the nonequihbrium DOS with increasing bias voltage 
[see Fig. 2(b)]. Furthermore, the nonlinear TMR displays 
a deep dip in the linear regime and changes its sign at a 
certain bias voltage. This large change in the TMR vs. 
bias voltage reflects the different behaviors of the Kondo 
resonance in the P and the AP configurations. 

IV. SUMMARY 

We have investigated the low-temperature, nonequi- 
hbrium properties of a spin-valve system consisting of a 
QD connected to two ferromagnetic leads in the Kondo 
regime. Based on the NCA approach we find markedly 
different behaviors in the equilibrium DOS when chang- 
ing the relative orientation of spin polarization. In the 
AP configuration, we find that a single Kondo peak al- 
ways appears through the whole range of polarization 
< p < 1, as just in a QD connected to two normal 
leads. Increasing polarization p can slightly suppress 
the amplitude of the peak. In the P configuration, the 
Kondo peak descends greatly and splits completely into 
two peaks even for a weak polarization as low as p = 0.2, 
leading to a steep valley with nearly a zero bottom. In 
both configurations the chemical-potential difference (the 
bias voltage) appears in the DOS via the the splitting of 
the Kondo peak into two peaks. Thus four peaks can be 



found in a moderate bias voltage for the P alignment. Of 
course the amplitudes of these peaks are suppressed by 
increasing temperature. 

Experimentally, we predict, based on the NCA investi- 
gation, that the different Kondo effects can be observed 
in transport through a QD by either linear or nonlinear 
measurements. For the AP configuration, the calcula- 
tions exhibit the usual temperature dependence of the 
linear conductance and a zero-bias maximum in the non- 
linear conductance, which are the conventional proper- 
ties of the Kondo-dominated transport through a QD. 
For the P configuration, however, we find a "hump" 
in the temperature-dependent linear conductance and a 
nonzero-bias maximum in the differential conductance. 
These peculiar behaviors are associated with the fact that 
the Kondo peak of the QD is split and shifted away from 
the Fermi level in the case of the P alignment. Further- 
more, we suggest that the TMR is a more effective tool 
to explore the different features of the Kondo resonance 
in different configurations. 

Finally, we point out again the applicability of the 
NCA approach to the FM/QD/FM systems. No spin 
splitting in the case of the AP alignment guarantees that 
the NCA is a reliable approximation for quantities in- 
volving the DOS. For the P configuration our numerical 
results for the DOS show satisfactory agreement with 
previous EOM predictions in the case of weak polariza- 
tion p = 0.2. Thus we believe that the self-consistent 
second-order perturbation approach provides some qual- 
itative features of the Kondo effect in the DOS as long as 
the polarizations of the leads are weak enough, which 
can serve to furnish a deeper understanding of trans- 
port properties in FM/QD/FM systems. Advanced NCA 
studies containing the vertex corrections are required to 
examine the validity of this approach. Work along this 
line is in progress. 
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FIG. 1: The total equihbrium DOS T{lo) in the AP con- 
figuration for different temperatures T/Fo — 0.01, 0.02, and 
0.04 and different polarizations p = 0, 0.2 and 0.4; and (b) the 
nonequihbrium DOS for T = 0.01 and p = 0.2. The QD in the 
FM/QD/FM system has a single bare-level energy = —4.0 
and an infinite on-site Coulomb interaction [/ — » oo. 
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FIG. 2: (a) The total equilibrium DOS in the P configuration 
for different temperatures T/To = 0.01, 0.02, 0.04, and 0.06; 
and (b) the nonequilibrium DOS for T = 0.01 and p = 0.2. 
The system is the same as described in Fig. 1. 




FIG. 3: The linear conductance G versus temperature at 
different polarizations p for the same system as described 
in Fig. 1. The thick linos correspond to the results for the 
P alignment, while the thin curve for the AP alignment at 
p = 0.2. Inset: The calculated TMR vs. temperature for 
p = 0.2 
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FIG. 4: The linear conductance G and the TMR versus the 
bare- level energy of the QD at T = 0.1. 




FIG. 5: The differential conductance dl/dV versus the bias 
voltage eV at different temperatures T = 0.02, 0.04, 0.06, 
and 0.1 in the P (a) and the AP (b) configurations, (c) The 
nonlinear TMR versus the bias voltage. 



